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Abstract. We present a concise introduction to a new spinwave expansion scheme for magneti- 
cally ordered Kondo lattice models. This is motivated by consideration of the ferromagnetically 
ordered phase of the "double exchange" system Lai_ v Ca x Mn03. A brief overview is given of the 
consequences of quantum and thermal fluctuation effects for the magnetic properties of the double 
exchange ferromagnet. 

INTRODUCTION 

The defining feature of transition metals is mixed valence — their ability to donate 
a different numbers of electrons to chemical bonds according to the detailed charge 
balance of their environment. In addition, magnetic transition metals like Manganese 
have a strong Hund's first rule coupling. This forces all electrons in the unfilled shell of 
the ion to align their spin, giving rise to a net local magnetic moment. 

Taken together, these two facts have important consequences for magnetism in doped 
transition metal oxides, where the transition metal ions can frequently be found in states 
with different valence and therefore different total spin. Where the electrons are also able 
to move from site to site, local moment behaviour normally associated with magnetic 
insulators can coexist with metallic conduction. In this article we will give a brief 
pedagogical overview of a systematic method of describing the magnetic and charge 
excitations of such as a system. The approach developed has the conceptual and technical 
advantage that it does not introduce any artificial distinction between the magnetism of 
itinerant and localized electrons 



FROM MATERIAL TO MODEL 

A good example of a mixed valent system with strong Hund's rule coupling is provided 
by the family of "colossal magnetoresistance" manganites Lai_ x Ca x Mn03. This com- 
pound has a perovskite structure, like the High Temperature superconductors, and the 
magnetic manganese atoms occupy sites on a cubic lattice — see Fig. U We can rewrite 
the chemical formula of this manganite as (La 3+ Mn 3+ )i_ x (Ca 2+ Mn 4+ ) x (0 2 ~)3. Strong 
Hund's rule coupling then means that for each unit cell the system contains 1 — x spin 
5 = 2 local moments (the Mn 3+ ions), and x spin S = 3/2 local moments (Mn 4+ ions). 
States with an integer number of electrons per Mn site cannot conduct, and the parent 
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FIGURE 1. Cubic lattice of Mn ions in Lai _ x Ca x Mn0 3 ■ 



compound LaMn03 is a spin 5 = 2 antiferromagnetic (AF) insulator. 

On doping with Ca, "holes" are introduced onto the Mn sites, so that the system now 
has a mixture of 5 = 2 Mn 3+ and 5 = 3/2 Mn 4+ ions. The crystal field at the Mn sites 
splits Mn J-electron levels into a ti g and an e g multiplet, as shown in Fig. ||. Where 
the e g electrons delocalize, this leads to a natural description of the material in terms of 
localized and itinerant electrons. 
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FIGURE 2. Crystal field splitting and Mn spin states in Lai_ v Ca T Mn03. 

The simplest model Hamiltonian we can use to describe such a system is the Kondo 
lattice model 

& = ~ f E c L c Ju -J^SiJi Si = ^cJjOapC/p (1) 

(ij)a i cc|3 

This model has already been much discussed during the course of this school, usually 
for small 5=1/2 and in the limit | / |<C t where J < is the AF Kondo coupling. Here 
we must consider the opposite limit of t <c| J | for ferromagnetic (FM) J > 0. We will 
consider only states in which the system is magnetically ordered, in which case it is 
reasonable to assume that the size of the local moment 5 ^> 1 . For simplicity, we neglect 
the degeneracy of the itinerant electron e g orbital. 




FIGURE 3. Paramagnetic state in the absence of inter site hopping (t = 0). 



Since t/J is a small parameter it makes sense to diagonalize the Hund's rule coupling 
term J first, and then to reintroduce the kinetic energy oc t. In the absence of any 
coupling between sites (t = 0), for any density of holes, the system must be an insulator 
with paramagnetic Curie law susceptibility — see Fig|3]. For finite t > the ground state, 
at least at large S, will be a metallic FM, since the band kinetic energy of the itinerant 
electrons is greatest when all the spins are aligned — see Fig This effect is was named 
"double exchange" by Zener [[]]], because in Lai_ x Ca x Mn03, itinerant electrons move 
between from one Mn site to another through connecting 02 p orbitals. Where such a 
FM groundstate is realized, it is therefore referred to as a double exchange ferromagnet 
(DEFM). 

In what follows we will develop a controlled expansion scheme for calculating the 
properties of the localized/itinerant electron system described by Eqn. |T]. Here we will 
only discuss the DEFM, but the methods developed are much more general and may be 
applied to any magnetically ordered Kondo lattice model. 
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FIGURE 4. A single double exchange bond in the limit / — > °°, S — > °°. 
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FIGURE 5. DEFM for t > and semiclassical mapping onto independent spin and charge excitations 
in limit J — > °°, S — > °°. 

SEMI-CLASSICAL TREATMENT 

The physics of the DEFM can be understood easily at a semi-classical level, by analogy 
with the problem of two local moments sharing a single electron. This was first solved 
by Anderson and Hasegawa O], who considered the Hamiltonian Eqn. [1] restricted to 
two sites. In the limit S — > °° we can define an angle between the two local moments 



(Fig. U>. The energy of the ground state for a given is given by 



£o = -y- ?cos Q) < 2 ) 

where an SU (2) rotation of the local spin quantization axis, together with the limit 
/ — > oo, has been used to project out those electron states which are not aligned with 
the local moment on each site. There are two contributions to this an energy, a term 
proportional to J which is the Hund's first rule binding energy of an electron on a single 
site, and a term proportional to t which is the kinetic energy gain which comes from 
sharing the electron between two sites. The factor of cos(0/2) comes from the SU(2) 
rotation. Using the relation cos0 = S1.S2/S 2 , and a standard trigonometric identity for 
double angles we can write the energy of this state as 

t - -> 

E » -— 2S1.S2 + const. (3) 

The kinetic energy of the itinerant electron generates an effective FM Heisenberg ex- 
change interaction Joe = t/AS between the two sites. The generalization of this semi- 
classical analysis to the lattice 0], gives an effective Hamiltonian for the DEFM 

Xeff = ^p+^a (4) 
with independent spin and charge excitations 

= ~t E ///; = ~ j de £ Sjj (5) 

<v> <<7> 

Here the sum runs over nearest neighbour sites, and fj is the creation operator 
for an electron aligned with the local moment on site i. The charge part of the Hamil- 
tonian describes electrons with a tight binding band dispersion £^ = —ztjk and the 
spin Hamiltonian has spinwave excitations with dispersion (x) q = zJdeS[\ — J q ] where 
Jk = (cos ft* + cosky + cosfc z )/3 is the structure factor for the cubic lattice of Mn sites 
in units where the lattice constant a = 1 . The size of the effective exchange interaction 
Jde is determined by the ground state average of the itinerant electron electron kinetic 
energy in the FM ground state of the model, and so is also proportional to t 

J oe = ^L(f!fj)o (6) 

(ij) 

In fact the energy scales of both spin and charge excitations must be be proportional to 
t, since the role of hopping is to lift the degeneracy between the many different charge 
configurations for a given doping and spin state. 

Using this effective Hamiltonian, we can understand both the metallic conduction of 
the FM phase of the CMR Manganites, and the fact that their spinwave spectrum, when 
measured in inelastic neutron scattering experiments, looks a lot like that of a FM near- 
est neighbour Heisenberg model. However we have neglected all quantum fluctuation 



effects, which means we cannot explain, for example, the large damping of spinwaves 
which is also observed in Neutron scattering experiments. The most important feature 
missing from this Hamiltonian is the interaction between spin and charge excitations, 
which is present when fluctuation effects are taken into account. It is this interaction 
vertex which we now calculate. However, in order to do this, we must first re-derive 
the semi-classical treatment of the DEFM in a rigorously controlled manner. Since the 
literature on the DE model is too large to review in this format, we will not make any 
attempt to discuss the alternative approaches which might be taken to this problem. 

NEW SPINWAVE EXPANSION 

Since the Hund's rule coupling / is much larger than the kinetic energy scale t, we begin 
by setting t = and considering the atomic limit of the Kondo lattice Hamiltonian Eqn. |T|. 
It is our aim to first diagonalize the local Kondo interaction, before reintroducing the 
kinetic energy terms. Physically, it makes no sense to distinguish between the magnetic 
excitations of the local moment S and the spin of the itinerant electron s, so we consider 
instead the total spin T = S + s. In terms of this coordinate, the Hamiltonian Eqn. [T] for 
a single sight can be written 

= -l[T(T+l)-S(S+l)-s(s+l)} (7) 

Here S is the size of the local moment, s is the spin quantum number for the itinerant 
electron, and the total spin quantum number T can take on values T = S ± s. There 
are four possible states. The lowest energy state has the eigenvalue E = —75/2. It 
corresponds to one itinerant electron being perfectly aligned with the local moment to 
give a total spin T = S + 1/2. There are two states with eigenvalue E = and total spin 
T = S. In these the itinerant electron orbital is either empty or doubly occupied. There 
is also a high energy state with eigenvalue E = J(S+ l)/2, and total spin eigenvalue 
T = S— 1/2. This has one itinerant electron perfectly "anti-aligned" with the local 
moment. 

Now we introduce electron operators which create these states, when acting on a 
vacuum in which the local moment has some definite z-component S z . We introduce 
one Fermi operator {/,/'} = 1 to describe high spin states, and another Fermi operator 
{PiP } = 1 to describe low spin states. These are illustrated in cartoon form in Fig. 
Next we need to introduce a bosonic description of the spin algebra of the total spin T . 
We do this by generalizing the usual Holstein-Primakoff transformation used to describe 
rotations to spins of a fixed length [£|] to the case in which the length of the spin is an 
operator. 
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FIGURE 6. Cartoon illustration of new operator coordinates. 
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FIGURE 7. Interaction between spinwaves mediated by density fluctuations of the electron gas. 



where 



T = S+ J \ (10) 



The physical meaning of these spin operator is illustrated in Fig. Finally, we ex- 
pand the operators {f,p,a} in terms of the operators {w,<i,a}, where [a, a 1 ] = 1 is the 
Holstein-Primakoff Boson for the core spin S. We do this order by order in the small 
parameter 1 / y/S. We must make sure that all the correct canonical (anti-)commutation 
relations are conserved by this transformation at each order. At the end of all this hard 
work, we are left with an inverse transformation 



45 / JlS 



■I 



p 1- 



45 



+ 



fa 



+ 0(l/5 3 / 2 ) 



(ID 



which can be substituted into the kinetic energy term in Eqn. [T] to give the interaction 
between spin and charge excitations. In the limit J — > °° the high energy "p" states can 
be neglected entirely. In terms of the new coordinates, the effective Hamiltonian Eqn. |] 
reads 



#P = L e kfl fk Xo = Y, ro ?4^ 



(12) 



The spin and charge modes now interact according to 
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(13) 
(14) 



Here the charge density : /'/ : is measured relative to the average density of itinerant 
electrons, and corrections to the vertex at o(l/5 2 ) have been suppressed. The form of 
the vertex is physically correct — spin and charge quantum numbers are independently 
conserved (as they must be in the limit J — > °°), and the strength of their interaction 
is proportional to t. The consequence of this vertex is an effective interaction between 
charge and spin modes at o(l/5 2 ), shown in Fig[7|. These results were first introduced 
in [Q]. Further details of the transformation, including terms up o(l/5 2 ), are given in 
[§and[[7p. 



RESULTS FOR THE DOUBLE EXCHANGE FERROMAGNET 

We now briefly illustrate some of the consequences of the interaction between spin and 
charge modes for the physical properties of the DEFM. All examples given are for a 
system with a density of x = 0.3 holes per site, spin 5 = 3/2 and hopping integral 
t = 1. These values are chosen to correspond with the FM state of Lai_ x Ca x Mn03 
with maximumal T c . Since the Hamiltonian is quadratic in Fermi and Bose fields, 
and the commutation relatations of all operators are well defined, standard diagramatic 
techniques can be used to calculate the spin and charge response functions of the model. 

In Fig. H the semi-classical spinwave spectrum oo q of the DEFM is shown, together 
with the renormalized spectrum when the leading quantum corrections to the spinwave 
dispersion at are taken into account. The effect of the interaction is large, especially 
at the zone boundary. The spinwave over all scale of spinwave dispersion is reduced, 
and no longer has the simple form associated with a nearest neighbour Heisenberg 
FM. Physically, this is because quantum fluctuations about the FM groundstate generate 
effective non nearest neighbour exchange couplings. 

Spinwaves excitations are exact eigenstates of the Heisenberg FM, but are damped 
in the DEFM. This damping arises because spinwaves can decay into lower energy spin 
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FIGURE 8. Spinwave spectrum of the DEFM at zero temperature. Upper solid line — semi-classical 
result throughout the Brillouin zone. Lower dashed line — result when leading quantum corrections are 
taken into account. 



excitations, dressed with charge density fluctuations (electron-hole pairs). This damping 
is also large, as is illustrated in Fig. [5] 




FIGURE 9. Damping of spinwaves in the DEFM throughout the Brillouin zone. 

Temperature corrections to the spinwave dispersion in a Heisenberg FM are exactly 
proportional to the zero temperature dispersion, and scale as T 5 / 2 at low temperatures. 
Temperature corrections in the DEFM are substantially enhanced, compared with the 
semi-classically equivalent Heisenberg FM, and are nolonger simply proportional to the 
zero temperature dispersion. The comparison of the two models is illustrated in Fig. [TC. 
Detailed results will be presented in a forthcoming publication [[7p. 



CONCLUSIONS 

We have developed a controlled large S expansion of the Kondo Lattice model, based 
on local eigenstates of total spin. This approach can applied to any magnetically ordered 
state, but is particularly convenient when working in the "double exchange" limit J ^> 
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FIGURE 10. Log-log plot of finite temperature corrections to the spinwave spectrum of the DEFM, 
together with those of the semi-classically equivalent Heisenberg FM. Temperature corrections Aco(<7, T) 
are normalized to bare dispersion co(g), temperatures to the scale of spinwave dispersion JdeS. 



t. This approach has been applied to the double exchange ferromagnet, which has 
been shown to behave in a quite different way from the Heisenberg ferromagnet when 
quantum effects are correctly included. 
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